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This paper proposes the design and implementation of predictive controllers 
for synchronous reluctance motor drive systems to enhance their dynamic 
responses. The predictive speed and current controllers in this paper are 
designed in systematic procedures. The predictive speed controller is 
implemented by using Laguerre function procedure. The Laguerre function 
is used to simplify the algorithm and to minimize the execution time of the 
digital signal processor. For predictive current controller, a finite control set 
method improves the current tracking ability. The measured currents are 
used to predict the future phase-current based on the motor model. The 
optimal control inputs of both predictive controllers are determined by using 
a cost function minimization method. Experimental results show the 
proposed drive system provides a wide adjustable speed range, from 2 1/min 
to 1800 r/min. It has better performance than a proportional-integral (PI) 
controller including fast rise time, which is 0.9 second, small steady state 
error, which is 0.32 r/min, and small current ripples. A 32-bit floating-point 
digital signal processor, TMS-320-F-28335 DSP, is employed to implement 
the control algorithms. 
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1. INTRODUCTION 


Nowadays, the using of power electronics technology and control technology for AC motors has 


improved the performance of AC drive systems [1], [2]. AC motors are classified into three types: i) the 
induction motor, ii) the permanent magnet synchronous motor (PMSM), and iii) the synchronous reluctance 
motor (SynRM). The induction motor was commonly used in industry due to its low maintenance and simple 
structure. However, the PMSM has slighter size and higher performance than the induction motor does. The 
permanent magnet material used in PMSMs, however, is very expensive. To solve this problem, an emerging 
trend in the industry is to replace the PMSM with the SynRM, in which no permanent magnet material is 
required. 

The SynRM has become more and more widespread due to its simple structure and no 
permanent magnet material is required [3]—[5]. The SynRM drive system is also very attracted compared to 
induction motor drive systems because its control strategies are more straightforward [6]. Additionally, the 
SynRM is free from rotor loss and slip, which provides higher efficiency and easier control than the induction 
motor. 
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The controller is an essential part for a motor drive system. There are various types of controllers 
that usually applied for SynRM drives. The simple way is using proportional-integral (PI) controller which 
the implementation is relatively straightforward. However, the tuning of PI parameter is more challenging for 
wide speed range [7]. Some advanced control of SynRM drive system have been examined by some 
researchers. Shyu and Lai [8] proposed multisegment sliding mode control for SynRM. The multisegment 
sliding mode controller was implemented for speed controller and the current loop used hysteresis controller. 
The method showed robust speed response but the ripple is huge. Lin et al. [9] proposed adaptive 
backstepping control for SynRM to improve speed tracking response. However, the result showed that the 
performance was good in middle to high speed. Senjyu et al. [10] proposed PI-based speed controller for 
SynRM with extended Kalman filter to achieve high efficiency control against parameter variation. 
Nevertheless, the method was too complicated when implemented in digital signal processor. 
Mahfoud ef al. [11] proposed direct torque control strategy motor drives using model reference adaptive 
system. The method is used for sensorless speed control that is independent to the stator resistor. 
Accordingly, the dynamic response is quick. 

Among the advanced control technique, the predictive controller has been successfully applied to 
industrial applications [12]. Generally, a model-based predictive controller (MBPC) requires a precise 
mathematical model of the motor [13]—[15]. Then, the future control signal trajectory should be predefined, 
and the output variable's future behaviors should also be optimized. A cost function minimization technique 
is added to calculate the optimal control input. As a result, the space-vector pulse width modulation 
(SVPWM) technique is not required here. The applications of the MPC include power converters [16]-[19] 
and motor drives [20]-[22]. Many previous papers have investigated the MBPC in either the current-loop 
control or the speed-loop control [23] but not including both of them in SynRM drive system. 

In this paper, predictive controllers are proposed for SynRM drive in current loop and speed loop 
control system. A systematic design of MBPC improves the dynamics performances for SynRM. The control 
algorithm proposed in this paper is executed by using digital signal processor (DSP). Therefore, the hardware 
circuit is easy to design. Compared to previous research [1]-[23], the ideas in this paper, which include 
investigation of the speed-loop and current-loop predictive controllers for a SynRM drive systems, are 
original ideas. 


2. PREDICTIVE SPEED CONTROLLER DESIGN 

In this paper, the speed loop of the SynRM drive system uses model predictive speed controller in 
which the Laguerre function procedure is employed. The predictive speed controller is designed based on the 
mechanical model of SynRM. The detailed model predictive speed controller design is explained in the 
following subsection. 


2.1. Discrete model of SynRMs 

The MBPC is a discrete-time controller. Therefore, the uncontrolled plant should also be presented 
in the discrete form as well. The discrete transfer function of the mechanical model of a SynRM can be 
represented as (1)-(3) [8]. 


Gp(z) = Wrm(Z) = (1 = raya Kt/Jmo ya (1) 


ig(z) S+(Bmo/Jmo) Z-g 


According to (1): 


-Bm 
peeele: Tmo") (2) 
Bmo 
-Bmo 
g=elm ™ (3) 


where Jo is the inertia of the SynRM, Bro is the friction coefficient of the SynRM, Ky is the constant 
torque, T, is the sampling interval, and Z is the z-transformation. Then, by using the inverse z-transform of 
(1), it is not difficult to obtain as (4), 


@,,,(n+ 1) = ga,,,(n) + hi, (n) (4) 
where W,m(n + 1) is the predicted speed; w,(n) is the measured speed; and i,(n) is the g-axis current. 


Predictive controllers for synchronous reluctance motor drive systems (Muhammad Syahril Mubarok) 


100 im) ISSN: 2088-8694 


2.2. Augmented model of SynRMs 
According to (4) and taking into account the different operations on both sides [24], one can derive. 


A@m(N + 1) = GA@,m(n) + hAig(n) (5) 


Where Aw,,(n +1) is the difference of predicted speed; Aw,,(n) is the difference of present speed; and 
Ai, (n) is the difference of g-axis current. From (5) is expressed as (6). 


Wrm(N + 1) = gA@,m(n) + AAig(N) + @m(N) (6) 
By subtracting ,, (n) of (6) on both sides, is defined as (7), 

Com (1+ 1) = gGA@rm(n) + hAig(N) + Cam M) (7) 
wheree,, (n+1)ande,, | (n) are the speed errors, and Lo, (n) is expressed as (8). 

Corm() = Wrm(N) — Wrm(n) (8) 
Then the new state variable vector is selected as (9) [24]. 

X(n) = [40;m(2) A@-m(n — 1) Aig (n—1) Brgy! (n)]* (9) 


By combining (5) and (7), yields. 


AW ;m(n + 1) _ {9 AWym(N) +i 
sat) =| le Crop, (2) + [,]4i000) (10) 
A new augmented state variable is defined as (11). 
_ [4@n(n + 1) 
X(n+1)= Coden (+ 1) at) 


And a new matrix and vector are expressed as (12) and (13). 


0 
el ae (12) 
= A (13) 


Substituting (11)-(13) into (10), one can define (14) [24]. 

X(n+1) = AX(n) + BAi,(n) (14) 
According to (14), the output is defined as (15) and (16). 

y(n) = CX(N) + @rm"(M) = Orm(7) (15) 

c=([0 1] (16) 
2.3. Using Laguerre function for predictive controllers 

A long control horizon can cause high computational numbers of the control signal Ai,(n). In this 


paper, the control horizon is used to capture the future control trajectory. As a result, the Laguerre function is 
applied here. First, Ai, (7) is represented as (17) [24]. 


Aig(n + p|n) = L(p) n(m) (17) 
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Where p is the predictive step, n(n) is the coefficient of Laguerre, and L(p) is the Laguerre function vector. 
The predicted state variable vector is obtained as (18) and (19) [24]. 


X(n+p|n) = APX(k) + p(p)" n(n) (18) 

p(p)" = Yaa A? **BL(d) (19) 
Where ~(p)’ is the transition matrix. Then by replacing (18) into (15), is defined as (20). 

Orm(n + p|n) = CAPX(n) + COP)" 1 + m(n +P) (20) 


Then, the optimal actuation can then be determined by minimizing the cost function. The cost function is 
expressed as (21) [24]. 


J= TPP, X(nt pln) OX(n+ pln) + EPS" Aig (n+p)? Ms Aig (n +p) (21) 


Where N,, is the prediction horizon, M, is constant and Q = C TC is the weighting matrix. Taking the partial 
derivative of (21), we can obtain the minimum cost function is defined as (22). 


sae = ie? 2Ata(n) (Zp? OCP) QP)" + Mi) +7, 2X(n) (Zp2, 0) 24”) (22) 
The new variables in (22) can be defined as (23) and (24). 

T=, 9) Qp(p)" + M, (23) 

A= 21 P(p) QAP (24) 
From (22), we can set rT to be 0, the optimal cost function is defined as (25). 

Oa = —IT-*AX(n) = n(n) (25) 


Finally, by substituting (25) into (17), the optimal solution of control input is derived as (26) and (27). 
Ai,(n) = —L(0)I1-1AX(n) = —Ring X(N) (26) 
Rinpe = L(O)l-1A (27) 


Where Rinpc is the state feedback control gain. 

In the real world, the SynRM control system should provide reasonable physical constraints on the 
state variables Ai, (n) and i, (n) to prevent saturation or damage to the hardware. The analytical solutions are 
required to identify the active constraints, which determine the control input Aig(n) and ig(7). As result, the 


mathematical modifications are needed when the constraints are imposed in the system. The constraints of 
Ai, (n) and i, (n) are expressed as (28) [24]. 


L < max 28 
nin ne aie 


Aig” (is ae < [im | 


The MBPC for speed-loop is redesigned when the state variables in (28) satisfy the constraints using 
Hildreth’s quadratic programming procedure [24]. In this paper, the constraint of the Ai, (7) is checked first. 
After that, the ig(m) is checked later to ensure that it is not exceeding the limit of the q-axis current. The 
MBPC for the speed-loop scheme is shown in Figure 1. The computations of the matrix and vector are 
simplified to simple numeric computations, which can be easily executed by a DSP. 
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Figure 1. Predictive speed control scheme of the SynRM drive 


3. MODEL-BASED PREDICTIVE CURRENT CONTROLLER DESIGN 
In this paper, the focus is on a three phase, balanced, wye-connected SynRM. The three-phase stator 
voltage equation is expressed as (29) [25]. 


Vin = Yglm + Lg SB + Cm (29) 


Where m denotes the a-b-c-phase, v,, represents stator voltage generated by the inverter, 7, represents stator 
resistance, i, represents the stator current, L, represents the q-axis self-inductance, and e,, represents the 
extended back-EMF. The SynRM is driven by a three-phase voltage-source inverter (VSI) where the input 
power is dc-link voltage. The predicted stator current in the discrete-time form according to (29) is 
represented as (30) [25]. 


ie (m+ 1) = (1-22) ipa) + Com) ~ ea) (30) 


Lq 


Where the superscript p denotes the predicted value; i®,(n +1) is the predicted value of the stator current; 
Vm (1) is the stator voltage; and e,,(n) is the extended back-EMF. The time delay compensation is taken into 
account for the model-based predictive current controller design. Then according to (30), the predictive 
current can be expressed as (31) [26]. 


TsTs 


h(n +2) =(1-% ) (n+ 1) + 2 (p(n +1) — em (n+ 1)) (31) 


Where i?,(n + 2) is the predicted stator current; v,,(n + 1) is the predicted stator voltage; and e,,(n + 1) is 
the future extended back-EMF. The stator voltage is set to be constant because the sampling interval is very 
short. Therefore, the stator voltage at the (1 +1) and (n) sampling intervals are assumed to be the same. 
Moreover, in (31), the future extended back-EMF e,,(n + 1) cannot be measured. The extended back-EMF 
can be expressed as (32) [27]. 


€m(R) © Cm (= 1) = Ym (n= 1) — (=) ig (0) + Lin (n - 1) (32) 


According to (32), assumed e,,(n +1) = e,,(n) due to the short sampling interval. Then (31) can be 
rewritten as (33) [27]. 


TsTs 
Lq 


i(n+2) =(1 J+ + ©(y,(n) — em(n—1)) (33) 
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The minimum value of the cost function is obtained using a simple mechanism, it is expressed as (34) [25]. 
k = |(Gx(n) — in(n + 2))| (34) 


Where the symbol * denotes the reference command. In (34), the cost function is composed of the predicted 
current error resulting from the deviation between the current command and the predicted current. The 
predictive current is obtained from the d-g axis currents, which is transformed into the a-b-c axis currents. 
Then according to (34), an optimal voltage vector can be defined as (35) [25]. 


k(n) Iso_s, = min{ k(n) |sy,--- CM) Is} (35) 
The voltage vector obtained from (35) is implemented to manage the six switches of the inverter at the (7) 
sampling interval. Then, the model-based predictive current controller is responsible for determining the 
switching state of the inverter, as illustrated in Figure 2. 


Cost function 


Sint , minimization and 

Stator curren F . 

aoe selection of the 
prediction | i?(+2) 


Inverter 
next switching | SiSet .| 
state 
a |b {ec 


a. 74 : 
offic 
bo 
4 Pi ————_ 


Using inverter switching states 


: 
To determine phase voltage of d-q to a-b-c 
the motor ' transformation 
Speed-loop 
controller 


oF, (1) ,_(0) 


Figure 2. The MBPCC scheme for the SynRM drive 


4. IMPLEMENTATION 

In the SynRM drives system is composed of two main components, the software and the hardware. 
The software is executed by a 32-bit floating-point TMS-320-F-28335 DSP, which is fabricated by Texas 
instruments. This DSP is employed to run the MBPC. The sampling interval of the speed-loop is 1 ms and 
the current-loop is 100 us. 

The block diagram of predictive controllers for SynRM drives is shown in Figure 3(a). The control 
algorithms in the software programs include the model-based predictive speed- and current-loop controllers, 
the d-q axis to a-b-c axis transformation, and a cost function minimization. The control algorithm starts with 
the speed command w;,,(7) is an input command. Then, the MBPC executes the speed-loop to obtain the 
optimal value of the i,"(n) while the ig*(n) is set as a constant. This i,"(n) is used for the current command 
in the current-loop. The d-q to a-b-c transformation frame is implemented because the current-loop control 
algorithm uses the a-b-c frame. By measuring the stator currents at the (nm) sampling interval, we are able to 
determine the values of the stator currents at the (n+2) sampling interval. Finally, a cost function 
minimization technique is used to determine the best candidate for the switching function. As the DSP only 
execute simple addition, subtraction, multiplication, and division. The computations of the matrix and vector 
are converted into simple calculations before the program is implemented. 

Figure 3(b) illustrates the implementation of the SynRM drive system. The hardware consists of a 
SynRM, a three-phase VSI, some current and voltage sensing circuits, a reshaping circuit of encoder, and an 
interface circuit. The SynRM used in this paper is fabricated by Reliance Electric Company, type PS6H5012. 
The motor parameters are as follows: Lg= 67.2 mH, Lg= 148 mH, R,= 2.0 2, 4-pole, 3.4 A rated current, 
220 V rated voltage, 560 W rated output power, and 1800 r/min rated speed. 
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Figure 3. The implemented predictive controllers for SynRM drives: (a) block diagram and (b) hardware 


5. EXPERIMENTAL RESULTS 

There are several experimental results shown here to validate the proposed MBPC for SynRM drive 
systems. The detailed parameters are selected as follows: the d-axis current command ij = 0.5, the initial 
value of Laguerre function L(0) = 0.71, the prediction horizon N,, = 1, and the current limitations including 
ane = —0.5A, Ai’ = 0.54, i = —4A,i7* =4A. The prediction horizon is selected as 1 to 
simplify the required computation of the model-based predictive speed controller design. 

Figures 4(a) and 4(b) show the predictive controller responses when the motor running from low 
speed to the rated speed. The different speed command with step input is given to evaluate the transient 
response. As can be observe from those figures, the proposed drive system has an adjustable speed range 
from 2 r/min to 1800 r/min with satisfactory performance. 
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Figure 4. Transient responses at different speeds: (a) low speed and (b) high speed 
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Figure 5 show the results of the comparison of the proposed MBPC and the PI controller. These two 
controllers operate with and without load. Then, based on the pole assignment technique, the parameters of 
the PI current-loop and speed-loop controllers are determined. In this paper, the PI speed-loop controller is 
chosen as K, = 50 and K; = 2. On the other hand, the PI current-loop controller is chosen as K, = 50 and 
K; = 5. Figure 5 demonstrates the transient response with step command at 800 r/min. Compared to the PI 
controller, the MBPC has a faster transient response, including faster rise time and no overshoot before it 
converges to speed command. As well, it is important to note that the integral gain in the PI controller causes 
an overshoot in the speed response. 

The rise time of the proposed MBPC is 0.9 second, while the PI controller is 1.6 second. The steady 
state error of the proposed MBPC is 0.32 r/min, while the PI controller is 6.6 r/min. The load disturbances 
responses at 500 r/min are shown in Figure 6 A 1 N.m external load are added to the drive system while it is 
running at steady-state. Once again, the proposed controller provides better performance than the PI 
controller. The MBPC has a small speed dip and a quick recovery time. Figures 7(a) and 7(b) show the 
measured speed responses to the sinusoidal speed commands using two different controllers. As can be 
observed, the predictive controller has better tracking responses than the PI controller for both positive half 
cycle and negative half cycle. As shown in Figure 8, different controllers are used to track current in the a- 
phase circuit. The proposed MBPCC can track the current command very well. Furthermore, the MBPCC has 
better current tracking responses and smaller current ripple than the PI controller. 

Figures 9(a) and 9(b) show the measured current trajectory in the a — 6 frame when using the 
predictive controller and the PI controller. In Figure 9(a), the MBPC provides satisfactory current tracking 
with little current ripple. However, the PI controller has greater current ripple and has worse current tracking 
as shown in Figure 9(b). Figures 10(a) and 10(b) show the measured current errors in the d — q frame. In 
Figure 10(a), the predictive controller has little current error. However, the PI controller has greater current 
error in the d — q frame, which is shown in Figure 10(b). Actually, the current errors in the d-qg frame that is 
shown in Figures 10(a) and 10(b) have DC offset. The DC offset is caused by the A/D converter does not 
convert the exact value at zero points. 
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Figure 7. Speed responses of a sinusoidal commands: (a) positive half cycle and (b) negative half cycle 
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Figure 10. Current errors in the d — q frame: (a) predictive and (b) PI 


The design and implementation of a model-based predictive speed and current controller for a 
SynRM drive system is presented in this paper. Experimental results indicate that the proposed controllers are 
capable of achieving fast transient response and load disturbances response. In addition, the proposed 
controllers achieve wide adjustable speed from 2 r/min to 1800 r/min. In addition, the proposed drive system 
tracks the current command very well with small current errors. Although the analysis is complicated, the 
realization of the controllers is quite easy by using a DSP. The proposed predictive controllers can be 
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implemented in industrial application due to their satisfactory performance. However, the computational 
issue of the implemented predictive controllers for SynRM in this paper can be simplified in the future. 
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